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We present results for the equation of state for pure SU(3) gauge theory obtained with a renormalization-group 
(RG) improved action. The energy density and pressure are calculated on a 16"^ x 4 and a 32'^ x 8 lattice employing 
the integral method. Extrapolating the results to the continuum limit, we find the energy density and pressure 
to be in good agreement with those obtained with the standard plaquette action within the error of 3-4%. 



1. Introduction 

The study of thermodynamic properties of 
QCD is crucial for understanding the early Uni- 
verse and relativistic heavy-ion collisions |l[]. The 
data is encapsulated in the equation of state 
(EOS). Recently the Bielefeld group determined 
the EOS in the continuum limit for pure gauge 
theory using the standard plaquette action|Q. 
Extending this result to full QCD will require 
the use of improved actions to compensate the in- 
creased computer power necessary for full QCD 
simulations. As a first step of such a program, 
we have studied the EOS for pure gauge theory 
with a renormalization-group (RG) improved ac- 
tion A summary of results Q is presented in 
this article. 

2. Simulation parameters 

The RG-improved action we use has the form[^ 
Sg = coJ2i^ ^ Hoop) + ci^(l X 2 loop) with 
Co = 1 — 8ci and ci = —0.331. 

We perform simulations on 16'^ x 4 and 32^^ x 8 
lattices, and also on symmetric 16"* and 32"* lat- 
tices, at about 10 values of /3 = in the 
range T/T^ « 0.9-3.5. We generate 20 000 to 
36 000 iterations after thermalization on asym- 



metric lattices, and about 10 000 iterations on 
symmetric lattices, where one iteration consists 
of one pseudo-heat-bath sweep followed by four 
over-relaxation sweeps. Errors are determined by 
the jack-knife method. 

3. Temperature scale and critical temper- 
ature 

We fix the temperature scale using the string 
tension of the static quark potential: ^ = 



(aV^)(/3c) 
aV^)(/3) 



*Talk presented by M. Okamoto 



For this purpose, we fit results for ay/a 
using a scaling ansatz proposed by AUtonl^, 

{aV^m = /(/?)( 1 + C2 a{/3f + C4 a{p)^)/co (1) 

where f{f3) is the two-loop scaling function of 
SU(3) gauge theory, and power corrections in 
the pseudo lattice spacing a(/3) = are in- 

troduced to incorporate deviations from two-loop 
scaling, with /?i an arbitrary reference point. 

In Fig. |l], we plot the critical tempera- 
ture in units of the string tension, Tc/^/a = 
l/{Nta^/a{(ic))^ for the RG-improved action to- 
gether with the result for the standard plaquette 
actionQ. Making a quadratic extrapolation in 
1/iVt, we find Tc/^ct = 0.650(5) in the contin- 
uum limit for the RG-improved action, which is 
3% higher than the value 0.630(5) for the stan- 
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Figure 1. Tcl\fo as a function of l/N"^. 

dard plaquette actionQ. The discrepancy may be 
caused by systematic uncertainties in the deter- 
mination of the string tension for the two actions, 
which differ in the details. 

4. Equation of state 

4.1. Integral method 

We calculate the energy density e and pressure 
p using the integral method B: 



P_ 

y4 



d/^'AS" 



3p 



rp4 



T^AS. 
dT 



(2) 



Here AS = N^{{S)t - {S)o), where {S)t and 
{S)o are the expectation values of the action den- 
sity S — Sg/N^Nt at finite and zero temperature, 
respectively. The beta function dfi/dT is deter- 
mined from the scale parametrized by (|l|). 

Our results for the pressure p are shown in Fig.^ 
together with those from the standard action. We 
find that the magnitude of the cut-off dependence 
for the RG-improved action is similar to that for 
the standard action, and opposite in sign. 

We extrapolate the energy density and pressure 
to the continuum limit, assuming a quadratic de- 
pendence on 1 /Nt as the RG-improved action has 
discretization errors of 0(a2). In Fig.| results of 
the extrapolation are plotted by solid lines, and 
are compared with those for the standard plaque- 
tte action (dot-dashed lines). As observed in this 
figure, results in the continuuum limit for the two 
actions are in good agreement with each other 
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Figure 2. Pressure for A^t = 4 and 8. 

within the estimated error of 3-4%. 

We note that the curves for the plaquette ac- 
tion in Fig.|| are obtained from a reanalysis of the 
data in Ref. Q| employing the same ansatz for 
the string tension (|l|) in order to avoid ambigui- 
ties arising from the choice of scale. In practice 
we found the changes in the pressure and energy 
density due to the choice of scale to be very small, 
being less than the statistical error of 1-3%. 

4.2. Operator method 

The energy density can also be calculated 
through the operator methodQ: 

- * [c,((5,>-(5)o)-c,((5,)-(5)o)](3) 



rp4 



5 



where Cs and ct are the asymmetry coefficients. 
The pressure can then be obtained with the sec- 
ond equation of (|2|). 

In Fig. |4| we compare results obtained with the 
integral and operator methods. The one-loop per- 
turbative values pO|] are used for the asymmetry 
coefficients. We observe that the two sets of re- 
sults are consistent with each other. The remain- 
ing discrepancy may well arise from the use of 
one-loop asymmetry coefficients; for the plaque- 
tte action, non-perturbative effects are known to 
be important 

In the high-temperature limit one may cal- 
culate the EOS by perturbation theory. The 
leading-order results for the EOS are shown by 
horizontal lines in Fig. ^. As has been noted some 
time ago[p^, the perturbative value for A'^t = 4 is 
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Figure 3. Equation of state in the continuum 
limit for the RG-improved action (sohd hues) and 
for the standard action (dash-dotted hues). 

very smah for the RG-improved action. Our nu- 
merical results are much larger than this value, 
at least in the range T/Tc < 3.5 explored in the 
present work. 

A possible source of the discrepancy is a break- 
down of perturbation theory due to the infrared 
divergence. Another possibility is that pressure 
and energy density decrease towards the pertur- 
bative values at high temperatures. This, how- 
ever, would require an unusual situation of a neg- 
ative AS' since the pressure is expressed as an 
integral of AS with the integral method. 

5. Conclusions 

Our continuum result for the EOS of pure 
SU(3) gauge theory obtained with an RG- 
improved gauge action shows a good agreement 
with that from the plaquette action. This pro- 
vides a concrete support for the expectation that 
continuum results are insensitive to the choice 
of lattice actions. We also found that the en- 
ergy density and pressure for Nt = A overshoot 
the perturbative high temperature limit. Under- 
standing the origin of this behavior is left for fu- 
ture investigations. 
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